Quantum Kinetic Theory of Collisionless Superfluid Internal Convection 
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Superfluids can transport heat via simultaneous opposite flows of their spatially interpenetrating 
condensate and thermal components. While this internal convection is usually described within 
Landau's phenomenological two fluid hydrodynamics, we apply quantum kinetic theory to a dilute 
Bose gas held beween thermal reservoirs at different temperatures, and show that the phenomenon 
also appears in collisionless kinetic regimes, and should be directly observable in currently feasible 
experiments on trapped ultracold vapors. 



Because the dynamical simplicity of dilute ultracold 
gases allows us to distinguish essential mechanisms from 
the complex details that abound elsewhere, they are 
ideal many-body systems in which to probe the still 
poorly understood foundations of non-equilibrium sta- 
tistical mechanics. Non-equilibrium time evolution in 
trapped quantum gases is a topic of high current interest 
[l|-l3|. Besides investigating non-equilibrium time evolu- 
tion, however, cold gas experiments can also impose a 
wide range of spatial inhomogeneities, offering the com- 
plementary opportunity to study quasi-stationary trans- 
port phenomena. Yet experimental samples are typically 
so small and dilute that, although the gas does equi- 
librate by scattering and exhibits nonlinear mean field 
dynamics, kinetics is 'collisionless,' in the sense that the 
mean free path of quasiparticles is larger than the sys- 
tem size. In the consequent absence of local equilibra- 
tion, most forms of transport cannot appear. But in this 
Letter we show that an unusual form of heat transport, 
known as internal convection^ should be observable even 
in a collisionless Bose-condensed gas. 

Internal convection has been observed in superfluid He- 
lium y, and is usually described within Landau's phe- 
nomenological two-fluid theory [sj : under a temperature 
gradient, the thermal and Bose-condensed fractions flow 
through each other, in opposite directions, with no net 
flux of mass. As if to rebut the 19th century conclusion 
that cold is only lack of heat, condensate flows as a flux 
of cold, from cooler regions to hotter. This striking yet 
simple form of transport is a convection in the sense that 
it involves motional separation of the fluid, but in Fourier 
rather than physical space. While identifying superfluid 
and normal components in Helium can be subtle, the rel- 
ative motion of condensate and thermal fractions in a 
dilute Bose gas is determined straightforwardly, accord- 
ing to the standard Penrose-Onsager definition [Q], from 
the single-particle reduced density matrix. Direct obser- 
vation of these two fractions in time of fiight expansion 
is a standard thermometry technique in experiments. 

Since internal convection does not involve spatial sep- 
aration, we will for simplicity consider a quasi-one- 
dimensional sample. Although one-dimensional many- 
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FIG. 1. (color online) Internal convection: A condensed 
Bose gas in contact with heat baths of different temperatures 
Tl > Tn. Condensate flows towards the hotter reservoir, 
non-condensate in the opposite direction. 



body physics can differ dramatically from that in higher 
dimensions, for finite samples there exists a (quasi-) 
condensate regime 0, 13] in which exactly similar mean 
field and quasiparticle theory governs both cases. This 
has been confirmed experimentally, for example in tests 
of Gross-Pitaevskii soliton dynamics in one-dimensional 
systems Q. Our results are therefore not special to one 
dimension, but generic. Moreover we derive them in the 
collisionless limit, but hydrodynamic internal convection 
is already well known. And we estimate readily observ- 
able flow speeds, approaching the range of mm/s under 
currently attainable conditions. Our conclusion is that 
internal convection is a robust phenomenon that should 
be observable even well beyond the realistic regime that 
we consider. 

Heat transport in a dilute Bose gas under a tempera- 
ture gradient has indeed already been analyzed quantita- 
tively in [lo| , though without explicitly identifying inter- 
nal convection. Since collisionless gas does not support 
local equilibrium with a smoothly varying temperature, 
however, we take a different approach. Rather than as- 
suming a temperature gradient, we let two reservoirs of 
thermal gas, with unequal temperatures T^, T^, interact 
with a Bose-condensed system cloud in different regions, 
as sketched in Fig. 1. Such a configuration could be 
achieved by trapping atoms of two different elements, or 
in a spinor condensate with spin-flipping collisions sup- 
pressed 



llj . One species of atoms serves as the system. 



and the other populates the reservoirs, with the three 
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samples spatially localized by species-dependent external 
potentials |12i] . Species conservation prevents particle ex- 
change between system and reservoirs, but collisions will 
still permit energy exchange 13[. Solving for the sys- 



tem steady state reveals non-equilibrium transport that 
includes internal convection. 

We identify this steady state by finding the time- 
independent solution to a quantum kinetic master equa- 
tion for the full density operator p of the system gas. 
We derive this master equation by tracing over the reser- 
voir sectors of the total Hilbert space, using the standard 
Born-Markov elimination method of quantum optics as 
applied to_ quantum gas kinetics by Gardiner and Zoller 



et al. 
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For technical convenience we consider 



a discretized model with periodic boundary conditions: 
a Bose-Hubbard ring of 2M sites, with the two reser- 
voirs coupled only to the diamatrically opposite sites 
771 = 0, M, applying the reservoir-plus-localized-mode 
theory of [19|, l20|. We then take a continuum limit of 
our model to describe continuous condensates. 

The Born-Markov master equation presumes that the 
local atomic destruction operator at each Bose- 
Hubbard site m has the interaction picture evolution 



(1) 



for some set of eigenfrequencies uj and operator- valued 
Fourier co-efficients dm(w). Assuming that the reservoir 
coherence time is long enough to resolve all the system 
frequencies w, we obtain [ilj 



p,H 



7 Mm{uj,uj')p 

m = 0,M 
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M„^{uj,uj')p = R{l3„i{uj - uj'))[pal^{uj')am{uj),dl^am] 
-R{l3m{Lo' - ll>)) [al^{uj')am{(^)p, aLoi 

1, z<0 

e'\ z > 



R{z) = 



where il is the system Hamiltonian, 7 is the effective 
scattering rate of the reservoirs (taken as the same for 
both of them, for simplicity), and /3o,Af — ^1 i}iBTh,R) 
the inverse temperature of the respective reservoir. The 
non-Hamiltonian terms in our master equation ([5]) pro- 
vide energy gain and loss, weighted by appropriate ther- 
mal factors in R{z), and with site-dependent phases con- 
tained in the Fourier co-efficients ao.M- 

For H we take a one-dimensional Bose-Hubbard model, 
which can represent bosons trapped in a lattice poten- 
tial or, in the continuum limit, a bulk sample. This 
discretization will also facilitate comparison with trans- 
port problems iii spin chains j22i] and strongly interacting 
quantum gases 
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J{atia„r+i + H.c. 



(3) 



Here J is the tunneling rate between neighboring lattice 
sites, n = N/{2M) is the average site occupation, and 
K/{An) is the on-site interaction energy. The system's 
chemical potential p may be tuned for convenience, since 
N = aJ^Om commutes with H. We will consider an 
eigenstate of iV, with eigenvalue N. 

The master equation ([2]) with Bose-Hubbard Hamilto- 
nian ([3|) admits no general analytic solution. But in the 
strongly (quasi-)condensed case where almost all of the 
N ^ 2M particles occupy a single mode, so that 1 
while K J, we can solve for the steady state p analyti- 
cally within a number-conserving Bogoliubov expansion. 
We begin with the Holstein-Primakoff transformation 



-h (4) 



which leaves A as a canonical lowering operator for 
N = A'' A, while the bk are canonical bosonic destruc- 
tion operators that all commute with A (and so with N). 
We then assume a total number eigenstate N N, and 
expand H in inverse powers of average site occupation n. 
With a subsequent Bogoliubov transformation. 




'}k = UkCk + VkCl^j^f^k 



(5) 



for canonical quasi-particle operators Ck , and coefficients 
Uk±Vk = (cjfe/rife)^^^^, H reads: 



j=0 
2M-1 



(6) 
(7) 
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Herein, the eigen-frequencies fifc 
tain the eigen-frequencies ujk 
interacting model. 

With this Bogoliubov expansion we solve for p = by 
perturbing in l/fi and 7/^2 1- Assuming that although 
the reservoir coherence time is long, it is not long enough 
to resolve the very small frequency corrections from Hi, 
we base our Born-Markov interaction picture on Hq, so 
that we have oj —>■ ±57^, (or 0) in ([1]) and with 



A 



V2MN 



Ck 



A 



d™(-l],) = .,e-M'_^4^_^ 



(8) 



d„(0) 



A 



1 - 



N 



Inserting this with (|4]) and ([5]) into ([2|), we find that our 
master equation describes the particle-conserving excita- 
tion and damping of Bogoliubov quasi-particle modes in 
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the Bose-Hubbard system, due to scattering of reservoir 
particles from sites to = and M . 

Expanding H, Aij^ and p in 1/Vn and in 7, 



3=0 
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(9) 
(10) 



leads to the following set of equations: 
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Equations (ITTT) - p^ include no terms of order or n^/^. 
Such terms are in principle possible, but vanish identi- 
cally, showing that although internal convection is a co- 
herent effect, it is not superradiant. 

As usual in a perturbative master equation, the ze- 
roth order equation (|lll) determines poo to be diagonal 
in the Hq basis, but does not fix its diagonal entries. 
Rather, the diagonal elements of the left side of 
vanish identically, determining poo through the diagonal 
terms of Mmo- This yields a canonical distribution for 
each mode, but with a temperature that depends on flk 
(albeit weakly for small temperature difference): 



Poo 



In 
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The solution to the next higher order equations can be 
expressed in the form pij = pooCij. The detailed eval- 
uation of the operator-valued factors is lengthy but 
straightforward and therefore shown elsewhere 2l| . Here 
it suffices to note that Coi and Cn include only combi- 
nations of odd numbers of quasi-particle creation and 
annihilation operators and Ck, whereas Ciq consists 
solely of even-numbered combinations. 

To interpret our perturbatively obtained stationary 
quantum state we compute the single particle re- 
duced density operator X, which acts in a single-particle 
Hilbert space: 



[XUi = ^Tr (palai) < m,l < 2M 



Again we expand in 7 and 1 
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(16) 



(17) 



Note that the Bogoliubov expansion of cim means that 
different orders of pij contribute to the same order in X. 



It is now straightforward to show that Xqo = 1/(2M) 
and Xqi = 0. To zeroth order in 7n, therefore, X has 
one constant eigenvector with eigenvalue of order 1 (the 
condensate), and many eigenvectors with eigenvalues of 
order 1/7V (the non-condensate). 

By introducing the compact notation 
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we can express the leading term involving the rate of heat 
exchange with the reservoirs as 
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= Tr[pn(/3^,, + A)] 



(19) 



There is no contribution from Tr/5i2 because by normal- 
izing Trpoo = 1 we have forced all higher pij to be trace- 
less.) 

Given the forms of the eigenvalues of X to zeroth order 
in 7fi, one can easily show that, to order 771, X^^f'' will 
perturb only the small eigenvalues of X and their eigen- 
vectors, while XiQ will alter only the condensate number 
and wave function by adding a homogeneous phase gra- 
dient representing a coherent condensate fiow. Although 
the explicit evaluation of A"io is lengthy, therefore, it per- 
mits unambiguous analytical computation of the particle 
and energy currents, with the contributions to each from 
condensate and non-condensate fractions. 

The particle current can be defined by recognizing the 
discretized version of the continuity equation: 
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(20) 
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We observe that the expectation value j„i = T^T^pjm is 
nothing but an off-diagonal element of the one-particle 
density operator X, and so although we have only evalu- 
ated X perturbatively, no further approximation is made 
in extracting the particle current j,n from it. 

Since we do not allow particle exchange with the reser- 
voirs, there are no sinks or sources of particle density, and 
so our time-independent p must have uniform jm = j, in- 
dependent of TO. And since our circular system's station- 
ary state also has reflection symmetry through the diame- 
ter, this J must vanish. By distinguishing the condensate 
and non-condensate fractions, however, we reveal inter- 
nal convection. The non-condensate current can then be 
computed by 



j^^ = 2JIm 



^10 



(m — 1) ,m 



(22) 



Apart from small distortions near to = 0, A/, j^^^ — I for 
< TO < M, and j^^ = -I for M < m< 2M. In other 
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words, there is a uniform current I of non-condensate 
particles from the hotter reservoir towards the cooler, 
around either side of the ring. From our explicit analyt- 
ical results for piQ we obtain 21 1 



M 



/ = ^Ex™-M,x™.(i-(-ir') 



(2M)2 



(23) 



where Xm,k oc cos are Fourier factors, and is the 
expected quasiparticle population of mode k. 

While XiQ^ is obtained from pio, which involves only 
the Bogoliubov Hamiltonian Hq, X^q is given by pn, 
which includes the post-Bogoliubov dynamics of Hi . Be- 
cause the total current j vanishes exactly, however, we 
know that the condensate current balances the non- 
condensate: 



£ -.^ 2 Jim 



(rn— 1) ,m 



(24) 



In this sense we have exhibited internal convection as a 
post-Bogoliubov effect, without needing to perform any 
post-Bogoliubov calculations (though these are nonethe- 
less performed in 2l|, and confirm the result). The only 
property of Hi that we have needed is its contribution to 
continuity. Furthermore, one can show from the fact that 
H is real that neither condensate nor non-condensate can 
have any stationary current of order 7° (2l| . The / given 
by (j23l) is thus the internal convection current. 

We can also derive / for the continuum limit, by in- 
troducing a lattice spacing a, and keeping the system's 
length 2Ma fixed as a — 0. For a finite /, we must 
also extend our model straightforwardly by allowing each 
reservoir to couple to L/a adjacent sites, for some fixed 
L less than the total system size, rather than to only one 
site. For small relative temperature difference AT/T, we 
obtain the continuum internal convection current [2l| 



7r5/2 AT fkuTV L 



120 T \ p J A/' 



-pc , 



(25) 



where p is the linear particle density, c the speed of sound 
in the gas, p oa the chemical potential, and A/ = 
{V^pay^ the reservoir particles' mean free path in the 
sample gas (assuming all atoms have equal mass). A 
numerical analysis shows that the discrete result (|23|) is 
already well approximated by the continuum limit (j25p 
for system sizes around 100 lattice sites. 

Eqn. (|25p is our main result. It includes no dependence 
on the total system size: ballistic transport in the colli- 
sionless regime means that the effect depends on the tem- 
perature difference, rather than the temperature gradi- 
ent, as it would in the hydrodynamic regime. Currently, 
ksT ^ p is not difficult to attain experimentally, so it 



seems possible to achieve an internal convection conden- 
sate flow velocity I/p not far below c, i.e. approaching 
the mm/s range, for sufficiently large condensates. The 
velocity of the much less dense normal component will 
be correspondingly greater, to obtain equal current. It 
should therefore be possible to observe internal convec- 
tion as separation of condensate and thermal fraction in 
time of flight expansion. 

L. G. is a recipient of a fellowship through funding of 
the Excellence Initiative (DFG/GSC 266). 



[1 
[2: 

[s: 
[4: 

[6: 

[7 

[8] 
[9] 

[10: 
[11 

[12 

[13 

[14 

[15 

[16 

[17 

[18 

[19 
[20 

[21 
[22; 

[23] 



lgilzQrhrk.u ni-kl.de| 

T. Kinoshita, T. Wenger, and D. Weiss, Nature 440, 900 
(2006). 

B. Fischer, T. Schumm, and J. Schmiedmayer, Nature 
449, 324 (2007). 

W. S. Bakr, A. Peng, M. E. Tai, R. Ma, J. Simon, J. I. 
Gillen, S. Foiling, L. Pollet, and M. Greiner, Science 329, 
547 (2010). 

Superfluidity and superconductivity, 2nd ed. (Hilger, Bris- 
tol, 1986) ISBN 0-85274-807-8. 

L. D. Landau and E. M. Lifschitz, Hydrodynarmk (Berlin: 
Akad.-Verl., 1991). 

O. Penrose and L. Onsager, |Phys. R ev. 104, 576 (1956). 

D. Petrov, G. Shlyapnikov, and J. Walraven, Phys. Rev. 
Lett. 85, 3745 (2000). 

W. Liniger, Physical Review(1963). 

C. Becker, S. Stellmer, P. Soltan-Panahi, S. Dorscher, 
M. Baumert, E.-M. Richter, J. Kronjager, K. Bongs, and 
K. Sengstock, Nature Physics(2008). 

E. Zaremba, T. Nikuni, and A. Grifhn, Journal of Low 
Temperature Physics 116, 277 (1999). 

H.-J. Miesner, D. M. Stamper-Kurn, J. Stenger, S. In- 
ouye, A. P. Chikkatur, and W. Ketterle, Phys. Rev. Lett.] 

82, 2228 (1999). 

L. J. LeBlanc and J. H. Thywissen, |Phys. Rev. A| 75, 
053612 (2007). 

J. R. Anglin and W. H. Zurek, |Pi^i. Rev. Lett.| 83, 1707 

(1999) . 

C. W. Gardiner and P. ZoUer, Physical Review A 55 

(1997) . 

D. Jaksch, C. W. Gardiner, and P. ZoUer, Physical Re- 
view A(1997). 

C. W. Gardiner and P. ZoUer, Physical Review A 58 

(1998) . 

D. Jaksch, C. W. Gardiner, K. M. Gheri, and P. ZoUer, 
Physical Review A 58 (1998). 

C. W. Gardiner and P. ZoUer, Physical Review A 61 

(2000) . 

J. Anglin, Physical Review Letters 79, 6 (1997). 

J. Ruostekoski and D. F. Walls, Physical Review A 58 

(1998). 

L. Gilz and J. R. Anglin, to be published. 

C. Mejia-Monasterio and H. Wichterich, The European 

Physical Journal-Special Topics 151, 113 (2007). 

T. Prosen and M. Znidaric, 

Journal of Statistical Mechanics: Theory and Experiment] 

2009, P02035 (2009). 



